Colloids electrophoresis for strong and weak ion diffusivity 
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We study the electrophoretic flow of suspensions of charged colloids with a mesoscopic method 
that allows to model generic experimental conditions. We show that for highly charged colloids their 
electrophoretic mobility increases significantly and displays a mobility maximum on increasing the 
colloidal charge for all salt concentrations. The electrophoretic mobility of highly charged colloids is 
also enhanced significantly when ion advection is dominant, leading to a strong heterogeneity in the 
local electrophoretic response especially at low salt concentration, when ion diffuse layers overlap. 
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Effective electrostatic interactions between charged 
colloidal particles in solution play a fundamental role de- 
termining the macroscopic phase and rheological prop- 
erties of colloidal suspensions that are pivotal for many 
applications in material science [H, In particular, the 
response of these systems to applied electric fields en- 
hances the degree of experimental control and tunabil- 
ity of suspensions down to the nanoscale and provides 
a natural means to design and operate nanodevices. As 
a result, electrokinctic techniques are pragmatically ex- 
ploited to control fluid flow at small scales, as for exam- 
ple in nano-fluidic molecular sorting 0,|3Li nano-medicine 
protocols i or micro-pattern assemSyl. 

A theoretical and quantitative knowledge of elec- 
trophoresis is still not complete as understanding the 
electrophoretic flow (EF) of the solvent and the charged 
components (ionized macromolecules, counter- and salt 
ions) requires the combined treatment of solvent flow cou- 
pled to the diffusion and advection of the dissolved elec- 
trolyte in the presence of charged colloids. The dynam- 
ical coupling between suspended charged components 
through hydrodynamics, leading to effective electrostatic 
interactions, does not allow for exact solutions; analytic 
predictions arc restricted to the linear regime and either 
numerical or simulation studies are required to assess the 
interplay between electrostatics and hydrodynamics that 
results in the EF. The relevant (length|time) scales in- 
volved in EF range from the nano- (meter | second), char- 
acteristic of the charge distribution around the particles, 
to the colloidal size itself, up to the micro- (meter | second). 
This disparity of scales complicates the computational 
modeling because molecular dynamics is practically lim- 
ited to the detailed resolution of the smallest, nano- 
scales. Mesoscopic modeling (MM) which includes 

the appropriate solvent dynamics at a coarse scale, over- 
comes the difficulties of handling multi-scale simulations 
of a charged colloidal suspension EF. MM has been re- 
cently used to study the non-equilibrium dynamics of a 
variety of complex charged heterogeneous systems [lol |. 
and therefore offers a valuable way to obtain theoretical 
predictions for a comprehensive EF, in primis the fun- 



damental EF of a charged spherical particle. 

In this Letter, we study the electrokinctic response of 
charged colloids and examine the relevance of non-linear 
coupling and ion specificity for different types of EF due 
to charge density deformations induced by applied elec- 
tric fields. We do so by taking advantage of a novel meso- 
scopic simulation technique that naturally includes ap- 
propriate boundary conditions and does not require any 
assumption beyond Poisson-Boltzmann. We then cap- 
ture the significant physical couplings between charges, 
colloids and solvent, accounting for the non-linear re- 
sponse of the electrolyte rearrangements to the applied 
field. Moreover, our model can be used under general 
experimental conditions of salt concentration, ion diffu- 
sivity and colloid surface charge density. 

A charged macromolecule in solution is surrounded by 
a counter- and salt ion cloud which screens the colloidal 
charge depending on the salt concentration, c^, and the 
solvent dielectric properties. Applying an electric field E 
causes the electrolyte and macromolecule to move. The 
limiting particle velocity at the steady state, vi , develops 
as a result of the balance between the electrostatic and 
viscous forces. However, the deformed charge distribu- 
tion of CO- and counter-ions around the macromolecule, 
which gives rise to an electric double layer (EDL) around 
the object, is hard to derive analytically even in thermo- 
dynamic equilibrium and usually solutions are known for 
low electrostatic potentials at the colloid surface (referred 
to as the zeta potential, C)j when linear electrostatics 
holds [H, Moreover, non- linear EDL distortion due 
to the external electric field leads to further difficulties 
in its analytic understanding. 

O' Brien and White (OW) [l^l have integrated nu- 
merically the electrokinetic equations for the EF of an 
infinitely diluted spherical colloid in the linear regime. 
Their results indicate that the particle electrophoretic 
mobility, ^ = vi/E, depends on the colloid radius a, salt 
concentration Cs and dielectric properties of the solvent 
at temperature T through the Debye screening length. 
Ad = (STrZbCsZ^)^^/^ (for a symmetric electrolyte of va- 
lence Z-|_, here z+ = ±1, in a solvent with Bjerrun length 



Ih = e'^ / (AnekBT) with e the electron charge and 
ks the Boltzmann factor) and predicting either mo- 
bihty saturation or a maximum for high C values. Using 
different mesoscopic models, Kim efc al. [13[ were able to 
confirm OW predictions [l^ for low Cs, while Lobaskin 
et al. [3, [3 showed that the electrophoretic behaviour 
of salt-free systems can be systematically mapped to a 
corresponding low-salt suspension. Both cases, however, 
address the restricted case of small ^ and Cg, where OW 
predictions are expected to hold. 

Our mesoscopic representation of charged colloidal sus- 
pensions builds on a continuous description of the elec- 
trolyte, characterized in terms of the anion and cation 
local densities, p±. These densities evolve according to 
the electrokinetic equations, which read 



^ = - V • pfei7 + V • i5fc[V/9, + e/?z,pfeV(^], (1) 
ot 



dpv 
~dt 



r]'V'^v-\/p.id + P^ezkPk'^(p, (2) 



V'^ip = --y^ezkPk (3) 

fe 

where Dk, Zk,k ~ +, — arc the diffusivities and valences 
of positive and negative ions, p,v,pid and 77 correspond 
to the solvent density, velocity, ideal pressure and shear 
viscosity, Lp is the electrostatic potential and e the sol- 
vent permittivity. Eq. [T] expresses ion mass conservation 
as a result of diffusion and advection, while the solvent 
motion, Eq. [21 evolves according to the Navier-Stokes 
equation for a viscous fluid accelerated by electrostatic 
forces due to local charge density. Finally, the Poisson 
equation enforces the electrostatic coupling between the 
charged species and the macromolecules. 

The solvent motion emerges from a discrete lattice for- 
mulation of Boltzmann's kinetic equation Q coupled to a 
discrete solution of the convection-diffusion equation for 
the dissolved charged ion species [l^. Hence, we regard 
the counter- and salt ions as scalar fields within Poisson- 
Boltzmann (PB) level The colloidal macromolecules 
are individually resolved, embedded on the lattice and 
coupled to the fluid throu gh a ppropriate kinetic rules ap- 
plied on their boundaries 17[. The hydrodynamic forces 
exerted by the fluid on the suspended particles, together 
with the electrostatic and dispersion forces determine the 
motion of the macromolecules in the fluid [l8| . The finite 
resolution of the colloidal particles on a lattice requires 
a proper calibration to identify the effective size where 
C, is consistent with the colloidal charge . The elec- 
trostatic potential drop around the colloids emerges con- 
sistently as a result of the ionic dynamics coupled to the 
fluid flow without further assumptions. Although such an 
approach disregards ion correlations, it provides a general 
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FIG. 1. Dimensionless electrophoretic mobility for a spherical 
colloidal particle, jl, versus the dimensionless zeta potential, C- 
Results from simulations (filled symbols) and 0W[12| (empty 
symbols) for ka = 0.5(8.0), triangles down(up). 



framework to address electrokinetics at weak and strong 
couplings and identifies the relevant competing physical 
mechanisms in charged driven fluids. 

We start by analyzing the electrophoretic mobility, ^, 
of a crystal of spherical colloids of charge Ze and radius a 
as a function of Cg, Q and volume fraction, $ = 47ra^/3L^, 
where L stands for the system size. In order to compare 
with experimental and numerical results, we introduce 
the dimensionless mobility, p. = ^UiSLs.^ and zeta poten- 
tial, C = The units in the simulations ensure that 
the relevant scales keep the right ordering and are prop- 
erly resolved on the lattice 120|. We take na = 0.5,8.0, 
with K = I/A_D, representative of low to high salt concen- 
trations [12| , covering two orders of magnitude in molar- 
ity, ~ 10^^ — 10^^ M, when mapping our simulation sys - 
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tcms to polystyrene spheres with a radius of 17 nm 
The applied electric field is tuned to remain in the linear 
response regime, E <C C/^d- 

Fig. [1] shows the dimensionless mobility, p, against C. 
For small values, C ^ 3, we observe an excellent agree- 
ment with OW theoretical prediction [l^ , which holds in 
the linear regime for low values of C- Our data show that 
increasing ^ leads to an enhancement of fl with respect 
to OW. Such a deviation increases with ( and with the 
amount of charge in the system up to 50% for narrow dou- 
ble layers, na = 8.0. Although the theoretical predictions 
from OW display a maximum outside the linear regime 
and only for thin double layers (na ^ 2.75), the maxima 
observed in Fig. [T] appear at higher zeta potentials and 
develop for any na. Therefore, these results indicate that 
care has to be taken when using OW curves to compare 
with experimental results. 

The model described reproduces quantitatively exper- 
imental data for low salt-free colloidal suspensions. In 
Fig. [21 we show that agreement with experimental data 
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FIG. 2. Dimensionless electrophoretic mobility ft for a col- 
loidal particle with charge Z = 30 and radius a — 4.0 versus 
colloid volume fraction $ (triangle left). No salt is added to 
the system. Experimental (sphere and square symbols) and 
simulation results (empty triangles) from Lobaskin et aJ.[l^. 
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FIG. 3. Steady state colloidal velocity vi/vo versus Peclet flow 
number for na — 0.5, 1.0, 8.0. Inset: dimensionless mobility 
/i versus dimensionless zeta-potential, ^, for na — 0.5, Pe 
~ 10"^ (filled circles) and Pe ~ 0.3 (filled triangles). 



is obtained for latex particles suspended in water [IJ, [21 1 
for colloid volume fractions $ > 0.001. For very low 
the ionic contribution due to water dissociation is rel- 
evant [l^ and the mobility dependence on $ vanishes. 
Fig. [2] also displays simulation results where counteri- 
ons are resolved individually [l^ . The agreement shows 
that the electrophoretic response observed experimen- 
tally can be captured through the counterion density and 
that ionic finite-size electrostatic correlations are sub- 
dominant. 

Even if in many situations individual ion resolution is 
irrelevant, its diffusion may affect significantly the elec- 
trokinetic response of charged colloids. We have ana- 
lyzed the effect of ion diffusion for symmetric monovalent 
electrolytes, D = D±, and have quantified the impor- 
tance of electrolyte diffusion to flow advection in terms 
of the Peclet number, Pe = avo/D, where the charac- 
teristic fluid velocity is induced by the colloid and reads 
vq = eZE/67Ti]a. Fig. [3] shows the dimensionless steady 
state velocity vi/vq for a colloidal crystal at $ 10~^ 
as a function of Pe for different salt concentrations c^, 
Ka = 0.5, 1.0, 8.0, covering from narrow to wide EDL at 
a fixed colloidal charge density, ^ ~ 5.0, in the non-linear 
regime. The results show that for strong ion diffusiv- 
ity the mobility is not affected by advection, while at 
larger Pe a significant increase in the colloidal mobility, 
departing from the high diffusivity regime at Pe = 0, is 
observed. The inset of Fig. |3l where we display /i as a 
function of the colloidal C for regimes where ion diffusiv- 
ity is either subdominant or relevant, indicates that the 
competition between ion diffusion and advection by the 
incoming flow becomes quantitatively significant in the 
non-linear regime. The mobility enhancement reported 
at small electric fields arises from ion advection, there- 
fore differing from the mobility enhancement associates 



to charge stripping under strong applied fields [2^. Al- 
though we cannot discard that a maximum in fl develops 
asymptotically, for the range of ^ accessible experimen- 
tally, C< 10Q,[i2, we do not observe a saturation of fi. 
Results in Figs. [T][2] correspond to the small Pe ^ 10~^ 
regime to compare with theories and experiments where 
diffusivity is dominant; however the results of Fig. |3] in- 
dicate that diffusion can become a relevant variable be- 
cause a wide range of Pe can be reached experimentally. 
For example, for a colloidal particle of radius a 0.1 
/im, charge density a 8 ■ 10~^C/to^ and an external 
field Eexp lO^V/m Pe ~ 10'^ for an elec- 
trolyte composed of small ions with typical diffusivity, 
D ^ 10~^m^/s, Pe can nevertheless increase one or two 
orders of magnitude imposing strong electric fields or in 
the presence of bigger, nano-sized ions j23| which induce a 
greater EDL disruption. To characterize the EDL degree 
of distortion as a function of Pe, we compute the dimen- 
sionless eigenvalues A of the charge density inertia tensor 
around the colloid center, = Ai/(Za^), i = x,y,z. If 
we apply the external electric field along the x axis, the 
axisymmetric structure of the EDL gives an asymmetry 
of the transverse eigenvalue as Pe increases. Table U 
reports the magnitude of the three eigenvalues for differ- 
ent values of Pe and c^. For all salt concentrations Ai 
increase with Pe, showing a larger spread of the EDL 



TABLE I. Charge intertia matrix eigenvalues. 
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portancc of ion difFusivity identifying an experimental 
regime of enhanced mobility and increased heterogene- 
ity of local electrophoretic response when the external 
driving dominates over ion diffusion. At small ion dif- 
fusivities, the overlap of EDL leads to a strong dynamic 
coupling between colloids, causing strong heterogeneities 
in the colloidal response. We acknowledge J. Faraudo for 
fruitful discussions, the lEF Marie Curie scheme (GG) 
and the Direccion General de Investigacion (Spain) and 
DURSI project (IP) for financial support under projects 
FIS 2008-04386 and 2009SGR-634, respectively. 
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FIG. 4. Dimensionless mobility fl versus volume fraction $ for 
a suspension of particles with overlapping EDLs, Pe ~ 1, 10^. 
Inset: probability distribution of particle velocities vi when 
$ = 0.32, na = 1.0 (overlapping EDLs). 



that results in a lower screening of the colloid charge, 
therefore raising its mobility at steady state. The dif- 
ferent values of Ai show that the axisymmetric charge 
distribution around the colloid develops a significant de- 
parture from isotropy only when ion advection drag be- 
comes dominant, at experimentally high Pe. 

Finally, we study colloidal suspensions at low salt con- 
centrations and at volume fractions when EDLs over- 
lap. We compare our results for C = 2.0 with analyti- 
cal predictions derived by Ohshima using single-particle 
charge distributions and linear electrostatics l2J] . Fig. |3] 
shows that deviations from Ohshima's theory develop for 
$ ~ O.I, in agreement with the theoretical expectation 
that EDL overlap for $ > Ka/{1 + kq)^. We unveil also 
the fundamental role played by advection-diffusion com- 
petition in EF, as Pe leads to an increase in the elec- 
trophoretic mobility. Moreover, ion diffusivity is seen to 
alter conjointly the dispersion of vi in the linear regime 
for significant EDL overlap, as depicted in the inset of 
Fig. m Understanding the physical mechanisms underly- 
ing the increase in with $ for highly diluted suspensions 

^, requires further 



with overlapping EDLs j2Ij. j> < 10 
investigation beyond current computational capabilities. 
The probability distribution function (pdf) of colloid ve- 
locities for $ ~ 0.32, Ka = 1.0 shows that on increasing 
Pe not only the mean vi shifts to larger values but also 
the pdf widens. As a result, an increase in ion drag due 
to EF leads to larger spread in local vi, enhancing the 
colloidal local electrophoretic response heterogeneity. 

In summary, we have comprehensively studied the elec- 
trokinetic response of suspensions of charged colloids, 
opening up new possibilities to fine tune and extend cur- 
rent electrokinetic techniques. We have shown that the 
colloidal mobility develops a maximum due to the com- 
petition between EDL distortion and the EF induced 
by the applied electric field. We have revealed the im- 



* 'giupponi@f fn.ub.esl 

^ ipagonabarr aga@ub.edu| 

[1] W. B. Russel, D. A. Saville, and W. R. Schowalter, Col- 
loidal dispersions (Cambridge University Press, 1991). 

[2] R. G. Larson, The structure and rheology of complex flu- 
ids (Oxford University Press, 1999). 

[3] W. Sparreboom, A. Van Den Berg, and J. C. T. Eijkel, 
Nature Nanotech. 4, 713 (2009). 

[4] R. Austin, Nature Nanotech. 2, 79 (2007). 

[5] F. Andre and L. M. Mir, Gene therapy 11, S33 (2004). 

[6] R. C. Hayward, D. A. Saville, and I. A. Aksay, Nature 
404, 56 (2000). 

[7] P. J. Hoogerbrugge and J. Koelman, Europhys. Lett. 19, 
155 (1992). 

[8] R. Benzi, S. Succi, and M. Vergassola, Phys. Rep. 222, 
145 (1992). 

[9] A. Malevanets and R. Kapral, J. Chem. Phys. 110, 8605 
(1999). 

[10] I. Pagonabarraga, B. Rotenberg, and D. Frenkel, Phys. 

Chem. Chem. Phys. 12, 9566 (2010). 
[11] J. Barrat and J. Hansen, Basic concepts for simple and 

complex liquids (Cambridge University Press, 2003). 
[12] R. W. O'Brien and L. R. White, J. Chem. Soc, Faraday 

Transactions 2 74, 1607 (1978). 
[13] K. Kim, Y. Nakayama, and R. Yamamoto, Phys. Rev. 

Lett. 96, 208302 (2006). 
[14] V. Lobaskin et al., Phys. Rev 
[15] B. Diinweg et al., J. Phys.: 

(2008). 

[16] F. Capuani, I. Pagonabarraga, and D. Frenkel, J. Chem. 

Phys. 121, 973 (2004). 
[17] A. J. C. Ladd, J. Fluid Mech. 271, 285 (2006). 
[18] B. Rotenberg, I. Pagonabarraga, and D. Frenkel, Faraday 

Disc. 144, 223 (2010). 
[19] G. Giupponi and I. Pagonabarraga, Phil. Trans. R. Soc. 

A 369, 2546 (2011). 
[20] M 
[21] N. 

M 

T. 

[22] M 

[23' 
[24; 



Lett. 98, 176105 (2007). 
Cond. Matt. 20, 404214 



Gates et al, J. Phys: Condens. Matter 16 (2004). 
Garbow et al, J. Phys.: Cond. Matt. 16, 3835 (2004); 
Medebach and T. Palberg, ibid. 16, 5653 (2004); 
Palberg et al, ibid. 16, 4039 (2004). 
Fixman and S. Jagannathan, Macromolecules 16, 685 
(1983). 

A. Martin-Molina et al. Soft Matter 5, 329 (2009). 
H. Ohshima, Theory of colloid and interfacial electric 
phenomena (Academic Press, 2006). 



